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Abstract 

Elastic scattering resonances occurring in ultracold collisions of either bosonic or fermionic po- 
lar molecules are investigated. The Born-Oppenheimer adiabatic representation of the two-body 
dynamics provides both a qualitative classification scheme and a quantitative WKB quantization 
condition that predicts several sequences of resonant states. It is found that the near-threshold en- 
ergy dependence of ultracold collision cross sections varies significantly with the particle exchange 
symmetry, with bosonic systems showing much smoother energy variations than their fermionic 
counterparts. Resonant variations of the angular distributions in ultracold collisions are also de- 
scribed. 
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I. INTRODUCTION 

While numerous advances have been made in the study of dilute atomic gases at ultra- 
cold temperatures, very little is yet known about ultracold molecular gases, particularly for 
strongly polar molecules. The first experiments with ensembles of trapped polar molecules 



are only recently reported Molecular gases are more difficult to cool than their atomic 
counterparts, and the temperatures at which effects of quantum degeneracy are apparent 
are typically much lower. Improved cooling and trapping techniques for cold molecules have 
emerged, including the development of decelerators [2] and the formation of molecules from 
ultracold atoms via photoassociation [3]. But most trapped molecular ensembles are still 
relatively "hot" and rather far from the quantum degenerate regime. 

Many theoretical studies of molecular collisions at cold and ultracold temperatures are 
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One 



now emerging in order to spur and assist these ongoing experimental efforts 
of the most striking outcomes of recent studies is the identification of large elastic scattering 
resonances in zero-energy collisions of virtually all molecules aligned in an external electric 
field Q . These resonances can be "tuned" with the electric field strength, resulting in a novel 
kind of zero-energy spectroscopy [101], and they have also been predicted in calculations 



of trapped ensembles Uj. In an earlier work, we demonstrated the universality of these 



resonances for bosonic collision partners [12]. In this paper, we provide a number of details 
of the resonance classification, and extend the work by comparing and contrasting with 
similar resonances in collisions between fermionic polar molecules. 

A few general observations highlight the importance that resonance formation and de- 
cay will play in the dynamics of ultracold molecular gases: 1) near-threshold collisions of 
polar molecules are "resonance rich", with resonances of shape, Feshbach/Fano and mixed 
characteristics forming in the long-range dipole-dipole field, and with detailed characteris- 
tics sensitive to shorter ranged interactions; 2) the resonances occur systematically and can 
readily be predicted by simple semi-classical quantization formulas; 3) the character and 
frequency of resonances differ distinctly between bose and fermi collision partners; 4) the 
angular distribution of scattered products in near- zero energy collisions is strongly correlated 
with resonance formation and decay; and 5) the resonances are universal and, with minor 
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field variations, can be realized in virtually all strongly polar molecules. This ubiquity of 
threshold resonance phenomena implies that complex few and many-body processes in ul- 
tracold polar gases are likely to emerge as experiments progress to the quantum-degenerate 
regime. 



II. LONG-RANGE INTERACTION OF POLAR MOLECULES 

There are many molecules (such as alkali-halide salts) which, in their lowest energy states, 
exhibit strong polar characteristics. However, in the absence of external fields, as a conse- 
quence of parity conservation, even these strongly "polar" molecules have zero dipole mo- 
ment. Molecular interactions may polarize such collision partners, even in the absence of any 
external fields, though this generally produces a relatively benign and isotropic (1/r^) van 
der Waals attraction at large intermolecular separations. The stronger (1/r^) anisotropic 
interactions of interest here require the imposition of an external electric field, coupling 
opposite parity states, and producing a generally field- dependent dipole moment 

m - -(§ )^ (1) 

for all molecules. 

The magnitude and sign of the "induced" dipole moment, /i, therefore depends on specific 



details of the E{£) Stark shift of the state of interest for any particular molecule 13|. In 
general, for zero field, two molecular states of opposite parity are separated by an energy 
gap, so that E{£) varies quadratically for small fields. At larger fields. Stark maps tend 
to a linear regime, yielding a maximal induced dipole moment fl ^ d = qsd, where q is 
some effective charge, and s is a measure of the charge separation along the field axis, and 
typically is a few Bohr in size. The magnitude of the electric field required to maximize 
the dipole moment /i depends principally on the size of the zero-field energy gap; molecules 
with ground E states tend to have relatively large gaps and often require experimentally 
unrealistic field strengths, while those in 11 or A states can have quite small gaps (owing to 
the phenomenon of A-doubling), producing maximal dipole moments with modest fields. 
For ground vibrational states of S molecules fi{S) can be estimated as 
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where B is the rotational constant of the molecule and z is the field direction. In Tabled we 
provide typical values of the intrinsic dipole moment d and the rotational constant B for a 
variety of S-state molecules |l6| . 

TABLE I: Important dipolar and rotational constants of S 
molecules. 2.928(-05) stands for 2.928 x 10"°^. The dipole 
length, D is calculated in the assumption that the molecule 
is completely polarized. 
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The wide variety of polar molecules available for trapping (once adequate experimental 
techniques are developed) is one of several reasons for the intensity of recent interest in the 
field. 

At intermolecular separations much greater than the charge displacement, r ^ s, molec- 
ular interactions reduce to the familiar dipole-dipole form 



y(f) = ^[di-d2-3(f-di)(f-4) 



(3) 



which, for a pair of identical molecules in the same field-aligned state, yields an equation of 
relative motion 



V^ + ^(l-3(r.if) 



(4) 



2M 

where z is the field axis and M is the reduced mass. As noted above, assumes its maximal 
value of d at large fields, but reduces to zero as £^ — > 0. 
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Equation (jl]) is easily written in dimensionless form, independent of fi and M, 

V^(f) = ErelV^(f) . (5) 



-V + lfl-3(r-i)2 



2 ■ ' r3 

providing that all distances and energies are measured in dipole units (d.u.) 

D = a^M/h^ 

■ (6) 

Ed = nyiM^fx^) 

The emergence of these natural length and energy scales, previously noted by many au- 
thors is striking when one recognizes that for maximal /i, D ~ 10^ — 10^ times larger 
than typical molecular length scales, while Ed ~ 10^^ — 10~^^ times smaller than typical 
rotational level splittings. Dipole-dipole interactions dominate the long range part of the 
interaction in any polar system, suggesting that observable characteristics of polar gases 
at sufficiently cold temperatures will be universal. Two types of long-range dipole-dipole 
effects can accordingly be identified: those that do not depend on the short-range bound- 
ary conditions and those that result from an interplay between short-range and long-range 
physics. In this work we model the short-range interaction with a spherically-symmetric 
"hard wall" boundary condition at r = r^ = Rq/ D. 



III. ADIABATIC REPRESENTATION 



Expanding the solutions of Schrodinger's Eq. (5) in terms of spherical harmonics seems 
a natural numerical approach with obvious advantages: the symmetries of the dipole-dipole 
interactions are clearly revealed in the spherical harmonic basis. The interaction conserves 
the parity of the state and the cylindrical symmetry allows for conservation of the quantum 
number m. The states of a given orbital quantum number / are coupled only to the two 
states of / ± 2, which leads to a useful tri-diagonal structure of the potential matrix. 

In a spherical harmonic representation, however, channel coupling is not small compared 
to the diagonal elements for the whole range of Z = . . . oo, and in the case of the s-wave, 
which is expected to contribute considerably at low energies, the coupling to the d-wave 
dominates the interaction. It is advantageous, therefore, to change the representation so 
that an efficient channel decoupling is achieved. For that purpose, we utilize the Born- 
Oppenheimer adiabatic representation, and diagonalize the angular part of the interaction 
together with the kinetic centrifugal terms. 
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For a given value of the interparticle distance r we calculate the eigenstates of the following 
angular operator Ha 

[Ha{r)]n> = l{l + l)5n' + \{lm\l-?,{f-zf\l'm) . (7) 

As the potential matrix elements 

= (/m|l -3(f ■i)2|/'m) = 



'1 _ 3 / {l-m)(l+m) I {l-m+l){l+m+l) \ \c 

y 21+ 1 V 21-1 2l+-i 

3 /((;+l)2-ni2)(a+2)2-m2) 



'2l+3\j" ' (21+1)121+5) ^l,l'-2 

/({i-l)2-m2)(/2-m2) r 
Y (2/-3)(2i+l) ^'.''+2 



(8) 



2i-l 

depend on m^, all the adiabatic states are doubly degenerate except m = 0, which is non- 
degenerate. 

Explicitly, we solve the following eigenvalue problem for each r 

Ha{r)(j)nif]r) = Vnir)(f)nif;r) . (9) 

As the off-diagonal part of the dipole-dipole interaction decreases at larger distances r — > cxo, 
all the adiabatic states 0„ approach corresponding spherical harmonics and the potential 
curves Vn/r'^ approach the centrifugal barriers /(/ + l)/r^. This allows us to classify the 
adiabatic states according to their long-distance behavior. For instance, in the bosonic case 
the lowest state approaches the s-wave, the next contributing adiabatic state has rf-wave 
asymptotic behavior. We can therefore label the adiabatic potentials vim{f)/r'^ and the 
corresponding adiabatic eigenstates with quantum numbers n ^ l,m according to their 
long-distance angular behavior. 

Constructing the adiabatic solutions of Eq. ([9]) by expansion in spherical harmonics, 

0^™(r;r)= Y: '<l>Lir)YUr) (10) 

F>|m| 

we are now ready to expand the solutions of Eq. in terms of the adiabatic states 

V^(f) = ^F;™(r)0;„(f;r) . (11) 

l,m 

(The prime on the sum in Eq.f fTOj) indicates that only the terms with / -|- /' =even are 
included.) The radial wavefunctions Fim{r) then satisfy the system of equations 

(-^^ + ^,vUr) - E)FUr) + Y.{<i'\r)^^ + P^;;^\r))F,m{r) = (12) 



TABLE II: Convergence of the elastic cross section with respect to the number of channels for 
the incident energy oi E = 5 x 10~^ d.u. . Results for adiabatic (ad) and partial wave (p.w.) 
representations are shown. The adiabatic representation demonstrates much faster convergence, 
especially in the vicinity of a resonance. 
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The long-range behavior of Fim{r) coincides with the usual partial wave amplitudes, since the 
adiabatic components (pimi^'-,"^) approach the spherical harmonics (Pijn{r) as r oo. 

We have solved the coupled-channel Schrodinger equation f[T^ numerically, propagating 
solutions from the hard-sphere radius to asymptotic distances r ~ 500D. We employ quintic 
splines and orthogonal collocations for discretizing the Schrodinger equation and match 
the numerical solution with asymptotic boundary conditions at the right end of the interval. 

Unlike the partial wave expansion, the coupling matrix elements p and q are not tridi- 
agonal in /. For low energy calculations, however, the adiabatic representation is still ad- 
vantageous, as one can see from Table [ITl The adiabatic representation clearly convergences 
rapidly, even in the vicinity of threshold resonances where the partial wave representation 
fails to converge with a given number of channels. 

Neglecting the off-diagonal channel coupling in Eq. (fT2l) yields an adiabatic approxima- 
tion, in which elementary one- dimensional potentials 

ViUr) = + pir\r) . (13) 

afford a simple, intuitive picture of the scattering. These effective potentials are shown in 
Fig. [1] for bosonic (a) and fermionic (b) molecules. All the potentials have a very strong 
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FIG. 1: (Color online) Adiabatic potential curves for bosonic (a) and fermionic (b) polar molecules. 
The curves are labeled by their asymptotic angular momentum H. and by their exact magnetic 
quantum number |m|. Notice a small potential barrier in the lowest fermionic adiabatic potential 

^1,0. 

short-range attractive well. The long-range behavior in all the channels is ^, excluding 
the lowest Vq^q- This lowest channel is attractive everywhere with a ^ attractive tail. All 
the other effective potentials are repulsive at large distances and have a single potential 
barrier in the transitional region between short-range attraction and long-range repulsion. 
The barrier heights provide important additional energy scales for the system. As long as 
the collision energy does not exceed the barrier in the second adiabatic channel, the lowest 
channel strongly dominates and, therefore determines the character of the dipole-dipole 
scattering. This serves to identify the "threshold regime" as corresponding to collision 
energies that are much smaller than the effective potential barriers. Another important 
feature of the adiabatic potential curves is the presence of a small, ~ 0.1 d.u., barrier in the 
lowest fermionic channel. This barrier results in substantial differences between fermionic 
and bosonic threshold scattering. 

Contributions of various channels to the total elastic cross-section for bosonic scattering 
(averaged over the electric-field direction), are shown in Fig. [2] a), for tq = 0.0108 dipole 
units. The uncoupled Vq^o and ¥2,1 channels dominate bosonic scattering, while Vi^o and 
Vi^i dominate fermionic scattering, for energies below ~ 130i?D, as indicated above. Note, 
however, that detailed convergence of the cross-section (even at zero energy) requires several 
adiabatic channels. Both the magnitude of the threshold cross section and the degree of 
convergence vary with the hard-sphere radius, as they are sensitive to resonance formation 
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in the inner wells of excited potential curves. This is illustrated in Fig. [3l where the elastic 
cross section (averaged over the field direction) is plotted versus hard-sphere radius at the 
near-threshold energy, E = 5 x 10~^ E^. Also note that the threshold cross-section varies 
by 4-orders of magnitude, depending on the specific details of the short-range scattering. 




u.uui ua;i u. I I lu UA( u.uui vy.ui u. i i lu luu 

C) E.E„ d) E,E„ 

FIG. 2: Energy dependence of the field-averaged scattering cross section, a) Bosonic off-resonance 
scattering, tq = 0.108. b) Fermionic off-resonance scattering, ro = 0.035. c) Bosonic resonant 
scattering, rg = 0.042873 d) Fermionic resonant scattering, rO = 0.05693 

The nature of these variations can be understood better after a closer look at Fig. [H Each 
of the adiabatic channels has a strong short-range attractive core which can support many 
bound states. When varying the short-range boundary conditions we effectively vary the 
number of bound states supported by each of the potentials. Each time one of the adiabatic 
potentials looses a threshold bound state we see a strong variation in the total cross section. 

Figure 2 therefore illustrates four distinct cases, depending on the presence of a near- 
threshold bound state and the system parity. The effect of a near-threshold resonance on 
the cross section energy dependence is different for bosons and fermions. Even though the 
threshold bound state strongly influences the magnitude of the cross section in the bosonic 
case, the energy dependence of the cross section remains qualitatively the same (Fig[2^,c). 
In the case of fermions, however, there is a striking difference between resonant and off- 
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resonant cases (Fig [2}D,d): the shape resonance formed in the lowest adiabatic potential 
is narrow enough to produce a rapid variation of the cross section. In the non-resonant 
fermionic case (Fig[2)D), when the collision energy approaches the barrier top of ~ O.ld.u. 
the scattering cross section also grows, although not producing as distinctive a feature as 
in the resonant scenario (Fig[2]i). While bosonic threshold bound states produce a cross 
section enhancement in a broad range of the cut-off radii tq, fermionic cross sections are 
strongly enhanced only in very close vicinity of the resonant conditions (i.e. the resonances 
are clearly narrower). 

The influence of the short-range cut-off radius on the resonances' formation can be un- 
derstood from simple semiclassical analysis, as discussed below. 




FIG. 3: Near-threshold cross section as a function of the cut-off radius for bosons and fermions. 
Positions of the different {l,m) resonances predicted from WKB match numerical results. 

IV. SEMICLASSICAL ANALYSIS 

One of the advantages of the adiabatic representation is that it provides a simple inter- 
pretation of the entire set of resonances shown in Fig. [31 
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Let us look at the case of bosons first (Fig. [8^1). Consider tlie positions of wide peaks 
corresponding to a near-threshold bound state forming in the lowest adiabatic channel. Such 
bound states exist if the Bohr-Sommerfeld quantization condition (at = 0) is satisfied 



Although the semiclassical condition is not strictly applicable to the calculation of the po- 
sition of the rightmost peak in Fig. [3^ (since the right turning point in the Vo,o channel 
effectively lies at r — 00), we have introduced a single fitting parameter 0o = O.OTSStt 
in order to reproduce the position of the rightmost peak correctly. Fitting the potential 
Vo,o as 2r^Vo,o ~ S^t an explicit expression for the values of tq that support a 

near-threshold bound state 



The universal dimensionless parameters a = 0.9586 and f3 = 0.265 are obtained from fitting 
the adiabatic potential curve, and Eq. f|T4l) then gives the positions of the peaks in Fig. [3^ 
up to 3 significant figures. 

Similarly, the semiclassical description is useful for understanding positions of other reso- 
nances that occur in higher adiabatic channels. All the resonances in the model system can 
be classified according to the asymptotic behaviour of the effective potential that produces 
the resonance, i.e. the corresponding (/, m) quantum numbers, and the number of bound 
states supported by the corresponding channel potential. If nonadiabatic coupling could 
be neglected, the resonances would have zero-width when approaching threshold because of 
the large dipole- modified centrifugal barrier in Fig. [H In fact, they dominantly decay into 
the lowest m = adiabatic channel through short-range nonadiabatic coupling and show 
up as Feshbach resonances in the total elastic cross section. Their positions can be found 
directly from the usual Bohr-Sommerfeld quantization procedure with no fitting parameters. 
Although it is more difficult to develop a simple fit to the higher adiabatic potentials, an 
expression similar to f|T^ provides an excellent fit to the positions of the other resonances 
obtained from analyzing the semiclassical picture numerically 

r'o'"'\n) = \ . (15) 

ao + ain + 02^"^ 

We have performed such fitting for the dominant channels contributing to the low-energy 
scattering and the results are presented in Table lllli 





4«2 



(14) 



(4a + P{7m + 4>o)){nn + 0o) 
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TABLE III: Fitting parameters for resonant short-range cut-off positions rg for the channels dom- 
inating at low energy. 



{l,m) 


ao ai 


02 


Bosons 


(0,0) 


0.2447 3.3817 


0.711496 


(2,0) 


8.17602 13.5762 0.622266 


(4,0) 


31.6553 21.2334 


0.63889 


Fermions 


(1,0) 


0.3635 3.567 


0.6844 


(3,0) 


10.873 11.809 


0.68644 


(5,0) 


32.579 20.671 


0.63113 



The case of fermions is very similar to the bosonic one, but the potential barrier in the 
lowest adiabatic channel makes the overall picture quite different: It makes the resonances 
formed in the lowest adiabatic channel comparatively narrow in both the energy and cut-off 
radius domains. 

V. SCATTERING ANISOTROPY AND THRESHOLD RESONANCES. 

Unlike central potentials, the threshold scattering of dipoles can be strongly anisotropic 
even in the case of identical bosons. Although the (0, 0) channel strongly dominates scat- 
tering at low energies, the degree of this domination depends on short-range boundary 
conditions. For instance, when the total cross-section is minimal, the (0, 0) channel does 
not contribute to scattering, and we can expect the cross section to be strongly anisotropic. 
This is illustrated in Fig. (jib)- At the minimum of the total cross section in the vicinity 
of the (2, 0) resonance (Fig. H^) the isotropic component of the cross section practically 
vanishes, and we observe a strong anisotropy both in total cross section i (as a function of 
the incoming wave direction 6) and in differential cross section (as a function of the incoming 
wave direction and the scattering angle). In the resonant case, however, when the total cross 
section is maximal (Fig. Hli), both total and differential cross sections become essentially 
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isotropic with strong domination of the (0,0) channel (s-wave scattering). 

A similar situation occurs for fermions. The angular distribution (which is always 
anisotropic) is dominated by the rapidly varying (1, 0) adiabatic channel, and the rest of the 
channels that are not that sensitive to the short-range physics. As in the case of bosons, 
we observe a rapid variation of the angular distributions in the vicinity of a (3, 0) resonance 
(Fig. [5]). An interesting feature of the fermionic scattering, however, is the presence of a 
special direction 9 = tt/2, i.e. collisions perpendicular to the field. Since the contribution 
of the (1,0) component vanishes in this special direction, scattering perpendicular to the 
polarizing field is not sensitive to the details of the short-range interaction. 

VI. SOME OBSERVABILITY NOTES 

Sensitivity to tq indicates a corresponding sensitivity to short-range physics, and there- 
fore, for a given electric field, to the specific molecular species, and some molecules can 
be expected to be dominated by either resonant or non-resonant adiabatic channels. The 
dipole moment induced by an external electric field, however, varies with the field until the 
dipoles are completely polarized. This makes it possible to observe, in principle, the series 
of resonances which we described in previous section. 

Assuming the electric field is large enough to polarize the molecule, but small enough 
not to perturb the short-range wave function, the only part of the intermolecule interaction 
affected by the field would be the dipole-dipole interaction. Changing the induced dipole fi 
by tuning the electric field we can effectively manage the dipole scales. The elastic scattering 
cross section a will, therefore, scale with the field S as 

a = D{SYa{E/ED{S)) 

where D{£) and Ed^S) are given by Eq. [61 As aforementioned, the explicit dependence of 
the induced dipole moment on the field £ depends on the molecular state. 

As soon as the induced dipole becomes big enough for the dipole length to exceed the 
short-range scale about 7 times {Rq/D ^ 0.14), the first bound state is formed in the Vo,o 
potential and the scattering cross section peaks. This is the first peak in the series ([HI). 
This first peak would allow an experimentalist to identify the unknown empirical parameter 
Rq and estimate positions of the subsequent peaks of the series f[T^ from the condition 
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Tq (1) = 7^1) • The number of peaks that can be potentially observed for particular 
molecular species is limited by their intrinsic dipole moment: a.s fi < d, there is a maximal 
possible dipole scale Dmax, and, thus, the number of peaks rimax can be estimated from the 
condition r^^''^\nmax) > Tp^- 

Using the parameters for ^'''LiF molecules given in Table [T] as an example and deliberately 
choosing Rq = 80 short-range cutoff parameter, we show the direction-averaged 

total cross section as a function of the external polarizing field in Fig. [61 For weak fields the 
induced dipole moment is small and, thus, the dipole units are comparable to the short-range 
molecular scale. 

Another interesting and universal feature of the dipole-dipole scattering cross section as 
a function of electric field is the trend of going up with the field. This trend can be seen 
much more clearly in the case of fermions which exhibits much narrower resonant variations 
of the cross section. 

VII. SUMMARY 

An adiabatic representation is shown to be useful in numerical calculations of the ultra- 
low energy dipole-dipole scattering as well as for the classification of resonances that emerge 
in such systems. For both bose and fermi collision partners, the major resonant phenomena 
are determined by the m = channels, especially the lowest channel in the m = set. The 
resonant states formed result from the long-range part of the dipole-dipole interaction but 
are sensitive to interactions that take place at shorter scales. The contribution of short-range 
interaction can be modulated by tuning the applied field, making it possible to observe a 
series of resonances. The angular distributions of scattered particles, at ultracold energies, 
depend sensitively on the resonance positions, and display variable degrees of anisotropy. We 
have also identified a special scattering direction in the fermionic case that is not sensitive 
to the short-range physics of the scattered system. 
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FIG. 4: a) Total cross section for two bosons averaged over field direction in the vicinity of the (0, 0) 
(diamond) and (2, 0) (circles) resonances. b)-e) Total and differential cross section dependence on 
the angle 9 = cos~^ k ■ z between the incoming wave and the field direction for bosons. Both 
the total cross section and the differential cross section demonstrate strong dependency on the 
short-range physics. The dipoles are polarized along the z direction, xz defines the scattering 
plane. 
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FIG. 5: Total and differential cross section dependence on the angle = cos~ k ■ z between the 
incoming wave and the field direction for fermions. Both the total cross section and the differential 
cross section demonstrate strong dependency on the short-range physics. Scattering perpendicular 
to the field is not sensitive to the short-range conditions, however. 
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FIG. 6: Electric field dependence of the near-threshold ^'^LiF elastic scattering cross section, in 
both (a) absolute and (b) dipole units. The minima of the cross section in the upper figure rise 
rapidly with the field as a result of rapidly growing dipole length. The fermionic cross section 
exhibits much narrower resonant variations of the cross section while demonstrating the same £^ 
growing trend. 



19 



